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We present a first principles study of NMR and spin orbit effects in the unconventional super- 
conductor Sr2RuC>4. We have calculated the uniform magnetic susceptibility, which agrees rather 
well with the experiment in amplitude, but, as in an earlier model result we found the calculated 
hard axis to be z, opposite to the experiment. We have also calculated the Knight shifts and the 
NMR relaxation rates for all atoms, and again found an overall good agreement, but important de- 
viations from the experiment in same particular characteristic, such as the Knight shift anisotropy. 
Our results suggest that correlations in S^RuCU lead to underestimations of the orbital effects in 
density-functional based calculations. We also argue that the accepted "experimental" value for the 
relative contribution of orbital polarization in susceptibility, fO-15%, is also an underestimation. 
We discuss the puzzling invariance of the the O and Ru Knight shift across the superconducting 
transition for all directions of the applied field. We show that this fact cannot be explained by 
accidental cancellations or spin-flip scattering, as it happens in some elemental superconductors. 
We also point out that large contribution of the dipole and orbital hyperfine field into the Knight 
shifts in S^RuCU, combined with the possibility of an orbital-dependent superconductivity, calls 
for a revision of the standard theory of the Knight shift in the superconducting state. 

PACS numbers: 74.70.-b 76.60.-k 76.60.Es 74.70.Ad 74.25.Jb 



I. INTRODUCTION 

The layered perovskite Sr2Ru04 is often considered to 
be one of the rare cases of spin-triplet superconductors pi 
(see, however, Ref.H). While there is convincing evidence 
against a conventional s-wave stated the evidence for a 
spin-triplet p-wave superconductivity rests exclusively^ 
on the NMR experiments. Specifically, the main (and, 
arguably, the only) argument in favor of a (particular) 
p-wave state is the fact that the NMR Knight shift for a 
magnetic field parallel to the RuC>2 layers does not de- 
crease in the superconducting state. Indeed, in a singlet 
superconductor the Cooper pairs have no net spin and 
cannot produce any Knight shift, so for applied fields 
smaller, in energy units, than the superconducting gap, 
as the number of unpaired electrons at the Fermi level 
decreases so does the Knight shift. On the contrary, 
among the p— wave states allowed in a tetragonal sym- 
metry, there is one where the spins of the pairs lie in 
the ab plane and are capable of screening an external 
in-plane field. Indeed, a simple theory^ predicts that 
the Knight shift for such fields remains constant across 
the superconducting transition. The latter appears to be 
well established, being confirmed for the 17 shifts, for 
the 101 Ru shiftS, and for the spin susceptibility as probed 
by neutrons^. 

The superconducting state in question, sometimes 
called a chiral p-wave state, is described by a vector order 
parameter d oc {k x ± ik v )z, and is also compatible with 
the /iSR relaxation experiments^ which indicate appear- 
ance of spontaneous magnetic moments below T c that is 



usually interpreted in terms of a superconducting state 
with Cooper pairs having nonzero orbital moments. On 
the other hand, it is not readily compatible with multiple 
indications of line nodes of the order parameter—, since 
such lines are not required by symmetry in this chiral 
p-wave state and can appear only accidentally. While 
such an accident is always possible, no superconductors 
are known so far where the order parameter would have 
nodes not required by symmetry, for an obvious reason 
that such nodes strongly reduce the total pairing energy. 

The confidence in the chiral p-wave state has been fur- 
ther shaken by the fact that the Knight shift in Sr2Ru04 
remains constant not only for the in-plane fileds, but in 
fact for any direction of the applied field 9 ! According to 
the existing theory, no superconducting state allowed for 
tetragonal crystals can have such a property. The sug- 
gested explanation was that the direction of d in a mag- 
netic field of 0.02 T changes? 9 - creating a state with the 
symmetry d—k z x [a "rotated" state d—(k z +ik y )y: is not 
allowed in a tetragonal symmetry^ and can appear only 
as a second phase transition]. This explanation, however, 
is rather doubtful for the following reasons: (i) such a 
state would have an additional node line and therefore 
loose a good deal of pairing energy and (ii) although in 
this state the spins of the pairs lie in the yz plane, there 
is no y — z symmetry (as opposed to the xy plane) and 
it is not a priori clear whether the magnetic susceptibil- 
ity of the Cooper pair will be the same as for the normal 
electrons. Finally, (iii) spin-orbital part of the pairing in- 
teraction, which keeps the spins in the xy plane, despite 
z being the easy magnetization axisii, is rather strong in 
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FIG. 1: (color online) LDA (NMTO) band structure (thin 
line) and pure t2 g bands (thick line), obtained by integrating 
out all other states (see Table U}. The Fermi level is set to 
zero. 



this material, and the field of 0.02 T (1.1 fieV or 0.013 
K) seems to be way too small to overcome it. 

These considerations have spawn several alternative 
suggestions for the pairing symmetry in Si'2Ru04, such as 
a chiral d-wave stated or a mixture of nearly-degenerate 
planar p-wave states^. On the other hand, examples are 
known where the magnetic susceptibility and the Knight 
shift do not vanish at T — even in conventional super- 
conductors: mercury^ tinii or vanadium. 15 Such cases 
are traditionally attributed to either spin-flip scattering 
due to spin-orbit coupling on point defect s 16 i 17 or sam- 
ple boundaries 18 , or to accidental cancellations of the 
Fermi-contact and core polarization contributions due to 
peculiarities of the electronic structured 

It is at this point clear that first principles calcula- 
tions of the NMR Knight shift and the relaxation rate are 
highly desirable, in order to gain better understanding of 
the physics of the NMR in this compound. Similarly, to 
access the applicability of alternative interpretations it is 
important to understand better spin-orbit effects (note 
that selection between different p-wave states is entirely 
controlled by spin-orbit coupling; see Refs. [Hl2fj|). In this 
paper we present full-potential highly accurate calcula- 
tions of the spin-orbit effects in Sr2Ru04 within the local 
density approximation (LDA) with and without external 
magnetic field, as well as first principles calculations of 
the NMR Knight shifts and the NMR relaxation rates, 
and analyze the possible ramifications. In particular, we 
analyze the possibility of an accidental cancellations due 
to electronic structure peculiarities, the possible domi- 
nance of the van-Vleck term, and the effects of spin-orbit 
coupling and correlations beyond LDA. 



TABLE I: Hopping integrals t lmn in mRy. The connecting 
vector is T = lax + may + ncz, where a and c are the lattice 
constants; the hopping integrals are tabulated up to the fourth 
neighbors. 
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II. ELECTRONIC STRUCTURE 

Sr 2 RuC>4 has the perovskite crystal structure of 
the cuprate superconductor La2CuC>4 (space group 
I4/mmm)i2I The primitive cell is tetragonal body cen- 
tered (lattice constants a = b = 3.87A and c = 12.74A); 
Ru is located at (0,0,0) and Sr at (0,0,0.353), while the 
two non-equivalent oxygens are 01 at (1/2,0,0) and 02 
at (0,0,0.1615) respectively. 

The electronic structure of Sr2Ru04 has been already 
studied by several author o 22 i 23 ' 24 by means of density 
functional theory (DFT) in the local density approxi- 
mation (LDA) or in the generalized gradient approx- 
imation (GGA). In the present work we adopt differ- 
ent LDA approaches: the Linear Muffin Tin Orbitals 
(LMTO) method in the atomic spheres approximation 
(Stuttgart LMT047 code)^ the N-th order Muffin Tin 
Orbitals (NMTO)^ method and the NMTO-based down- 
folding approach (Stuttgart NMT047 code), which al- 
lows to calculate hopping integrals and Wannier func- 
tions, and the linear augmented plane wave (LAPW) 
method (WIEN2k code)^ 7 - Our results for the band 
structure (Fig. [1] and Fig. [5]) are consistent with pre- 
vious calculations. The NMTO bands (Fig. []} are ob- 
tained from the self-consistent LMTO potential, follow- 
ing a standard procedure^ 

The tetravalent Ru has 4 electrons in the d shell. The 
cubic crystal field splits the d levels into 3-fold degener- 
ate tig and 2-fold degenerate e g states; the bands at the 
Fermi level are thus ti g bands 2/3 filled. DFT calcula- 
tions show that these three ti g bands can be divided into 
a wider xy band, almost two dimensional, and two nar- 
rower almost one-dimensional xz and yz bands. In Fig.[T] 
we show the band structure obtained by using the NMTO 
method. The tig bands (thick line) were obtained by in- 
tegrating out all the degrees of freedom except for the ti g ; 
by means of this first-principles downfolding procedure 2 ^ 
we could construct Wannier functions and a real space 
Hamiltonian for these bands^. The corresponding hop- 
ping integrals are tabulated in Table U up to the fourth 
nearest neighbors; farther hoppings are tiny and can be 
neglected. The tables show that the inter-orbital hy- 
bridization is very small; thus the bare band dispersion 





FIG. 2: (color online) LDA (LAPW) band structure wil 
spin-orbit effects. The lines represent the actual band struc- 
ture (spin-orbit effects are small on this scale; the "error bars" 
indicate the size of the spin-orbit induced shift of each state. 
In the central panel we show the length of an "error bar" cor- 
responding to a shift of 500 K (w 43 meV). The horizontal 
line is the Fermi energy. 



can be written as 
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for each band. These bare bands are slightly modified by 
the hybridization terms that take form 



FIG. 3: (color online) LAPW Fermi surface, painted accord- 
ing to the spin-orbit splitting at each point (see the color bar 
for the scale in the temperature units). 



Since the bands at the Fermi level have mostly Ru ti g 
character with some admixture of O p, it is likely that 
NMR is also dominated by Ru d and O p contributions. 
In addition the polarization of the core electrons due to 
the d shell will be sizable for Ru. NMR results will be 
discussed in session IV. 



III. STONER FACTOR AND SPIN-ORBIT 
INTERACTION 

In the the spin density functional theory, the static spin 
susceptibility is formally exact, and can be calculated as 
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where t lmn are the hopping integrals tabulated in Table 
HJ This Hamiltonian yields a Fermi surface consisting of 
an electron-like sheet (7), which is neraly cylindrical in 
shape, and four crossing planes, reconnected into a hole- 
like (a) and an electron-like (/?) tetragonal prisms by the 
weak hybridization xz — yz. The Fermi surface obtained 
with the NMTO method is in good agreement with pre- 
viously reported results based on LAPW (Ref. 24). The 
LAPW Fermi surface is also displayed in Fig. O qualita- 
tively it is very similar to the results obtained in Ref. |24 
Differently from previous calculations here we include 
spin-orbit effects; these results will be discussed in Ses- 
sion V. 



Xs 



dM s (H)/dH, 



(2) 



where M S {H) is the spin magnetic moment induced by 
an external field H. This is related to the non-interacting 
spin susceptibility, which is simply the density of states 
(DOS) in appropriate units, xo = H%N(0), where is 
the Bohr magneton. The renormalization is quantified in 
terms of the Stoner factor, I : 



X,=Xo/0--IXo/2). 



(3) 



(in the literature one can find also an alternative defi- 
nition of the Stoner factor, per two spin instead of one, 
which differs from ours by the factor 1/2 in the denom- 
inator). One can also calculate the orbital part of the 
magnetic susceptibility as 



XL = dM L (H)/dH, 



(4) 



as long as the spin-orbit interaction is included in the 
calculations. 

The most straightforward and the most accurate way 
to compute I is by performing self-consistent LDA cal- 
culations in an external field. For this purpose we have 
selected the full-potential linear augmented plane wave 
method as implemented in the WIEN2k+lo package. 27 
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TABLE II: Static susceptibility in 10 4 emu/mol, as obtained 
from LAPW band structure. 
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This package allows for self-consistent calculations in an 
external field, including interaction of the field with the 
orbital moment. Two sets of calculation were performed, 
for a field along z or x axes. Particular care has be taken 
to provide high level of convergency: up to 729 inequiva- 
lent k-points (17x17x17 mesh), RKMAX up to 9, states 
up to 3.5 Ry included in the second- variational diago- 
nalization of the spin-orbit Hamiltonian. Local orbitals 
were included in the basis for both Ru and O, and a 
relativistic local p-orbital was added for Ru, to ensure 
the convergence of the spin-orbit calculations. The re- 
sults are shown in Fig. [4] and Table [TTJ The following 
observations are in place: (1) the spin susceptibility is 
anisotropic, with z being the hard axis, (2) the orbital 
susceptibility is nearly isotropic and (3) there is a meta- 
magnetic transition at a field between 75 and 100 T to a 
state with a magnetic moment that extrapolates to R^O.l 
/is/Ru in zero field. While this metamagnetism is obvi- 
ously inobservable, it is probably related to the metam- 
agnetism observed in some Sr 2 Ru04-based alloys. In this 
paper we will not discuss this metamagnetism and will 
concentrate on the susceptibility. The low-energy slopes 
of M(T) curves yield the results (in 10 -4 esu/mol) listed 
in Table HD 

We see that the calculations reproduce well the over- 
all scale of the susceptibility but predict an oppo- 
site anisotropy compared to the experiment, albeit of 
the same scale (6% vs. -11%). Note that similar 
problem was encountered by Ng and Sigrist in their 
model calculation^ They ascribed the discrepancy to 
anisotropy of the orbital susceptibility, which, they 
pointed out should be larger in the z direction because of 
the possibility to form an xz + iyz state out of degenerate 
xy and yz orbitals. Indeed, we find such an anisotropy 
in our calculated orbital susceptibility, but it is way too 
small. Yet we believe that Ng and Sigrist's conjecture was 
right and the reason we do not see this effect is due to the 
fact that the density functional calculations routinely un- 
derestimate orbital effects (overestimate orbital moment 
quenching). In the experimental literature one can find 
variety of estimates for the ratio Xorb/xtot, from 10%* 
to 17%22, while our result is 6%. Increasing our orbital 
susceptibility by a factor of two or three, however, still 
does not reproduce the experimental anisotropy. This 
strongly suggests, in accord with Ref. [2(| that the sep- 
aration of the orbital and spin parts in the experiment 
was inaccurate and in reality the orbital susceptibility 
is larger than usually quoted. Indeed, the separation of 
the orbital and spin susceptibilities in the experiment is 



based on the assumption that the orbital part is totally 
temperature independent, which is not true for metal- 
lic systems where a non-negligible orbital susceptibility 
comes from the states at the Fermi level (see, e.g, Ref. 
l3l| ). Note that the anisotropy of Xorb extracted in this 
way is 0.40xl0 -4 emu/mol^, 2.5 times smaller than the 
observed anisotropy of the total susceptibility 11 . We will 
return to this issue again later when discussing the NMR 
experiments. 

The calculation reported in this section yield a DOS 
of 50 states/Ry/cell. If interpreted in terms of Eq. ([3]), 
this results in the Stoner factor / = 0.46 eV, and the 
renormalization coefficient R = Xs/xo = 6.2. In the next 
Section we will be using LMTO calculation which give 
slightly smaller DOS, namely 45.6 states/Ry/cell. We 
will therefore use such R that renormalizes the LMTO 
DOS to the accurate LDA Xs, that is, R = 6.8. We 
point out that this value is about twice larger than the 
Stoner enhancing factor we obtained directly from self- 
consistent LMTO calculations in a magnetic field. This is 
not surprising; in our experience 32 allowing non-spherical 
variations of the spin density practically always results in 
more magnetic solutions. 



IV. NMR KNIGHT SHIFTS 

Next we calculate the NMR Knight shifts. These shifts 
originate from the hyperfme interaction — 7„I-H, between 
the nuclear magnetic moment — 7I and the hyperfme field 
H produced at the nucleus site by the electrons. The 
hyperfme field operator H is the sum of a core (H core ), 
a contact (H c ), a dipole-dipole (H d ) and an orbital (H°) 
term. In order to calculate the NMR shifts we used two 
different approaches. 

The first is described in Ref. |33j and based on the 
LMTO-ASA method. Here we just summarize the most 
important steps. The bare Knight shifts are calculated 
from the LDA band structure by using linear response 
theory and without introducing any external magnetic 
field. This means, e.g. for the spin terms, K = 2fi B Tr('\ 
\NH\ |), where we introduced the density-of-states op- 
erator (ks|AT|kV) = <5 ss /(5kk'(5(ek — Ep); here s and s' 
are spin indices, ek the band energy and Ep the Fermi 
level. This approach allows to clearly distinguish the 
contributions of the different orbitals and of the different 
interactions and thus makes the physical interpretation 
easier. The Stoner factor and the core shifts are obtained 
separately, i.e. from a self-consistent LSDA calculation 
with an external magnetic field, B ext . The latter pro- 
duces a core polarization m(r); the core shift is the ra- 
tio B core / B ext between the (core) contact hyperfme field, 
B C ore = /is(87r/3)m(0), and the external magnetic field. 
About 450 k points in the irreducible Brillouin zone were 
needed in order to obtain well converged results. We 
tested several different LMTO set up^i, with consistent 
results. The different contribution of conduction elec- 
trons to the Knight shift are displayed in Table IIII1 
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TABLE III: Knight shift, K a in %. The label a =x,y,z indicates the direction of the external magnetic field. The first panel 
presents the results obtained with the LMTO method in the scalar relativistic approximation; the Stoner enhancement factor 
R = 6.8 is used for the renormalized values. The second and the third panels are the results obtained with LAPW without 
and with spin-orbit (SO) interaction; the Stoner factor is already included by construction. The last panel lists available 
experimental data. For Ol the shifts for the (1/2,0,0) site are presented. 
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The second approach is the one implemented in the 
LAPW Wien2k packaged All the contributions to the 
Knight were obtained from the ratio between the hyper- 
fine field B and the external magnetic field B ext from 
which B originates. Clearly, all the Knight shifts calcu- 
lated in this way are already Stoner renormalized. The 
Knight shifts were obtained from the highly accurate 
LAPW band structure, and the effects of spin-orbit could 
be accurately analyzed. The results are again displayed 
in Table [ml 

Let's examine first the results obtained for the Ru 
isotropic Knight shift, K iso = (K x + K y + K z )/3 . We 
found that the s-electron contribution from the conduc- 
tion bands to Ki SO is negligible; this is understandable 
because N s , the Ru-s projected density of states (see 
Table |IVJ) is negligible at the Fermi level (N s /N t2g « 
2 x 1CP 4 ). The main contribution is the core polar- 
ization (the dipole-dipole term is sizeable, but purely 
anisotropic). 

At this point, it is interesting to compare the core- 
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FIG. 4: (color online) Orbital and spin magnetization 
(LAPW). 
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polarization Knight shifts in LMTO and LAPW (no spin- 
orbit), Table [Till As already explained, in both LMTO 
and LAPW the core polarization shifts were obtained 
from self-consistent electronic structure calculations in 
an external magnetic field, and are therefore renormal- 
ized by the Stoner factor. Still, the LMTO yields a core 
shift a factor 1/2 too small; this is consistent with our 
results on the Stoner enhancing factor R, which turned 
out to be a factor 2 smaller in LMTO than in LAPW. 
For this reason, in calculating the Stoner renormalized 
shifts for the LMTO (the rightmost column of Table ITTT]) 
we took into account this discrepancy and multiplied the 
core term by this missing factor 2. This is is an impor- 
tant technical point and we would like to reiterate it: 
In LAPW calculations all calculated contribution are by 
construction already Stoner- renormalized. In LMTO, the 
core polarization is renormalized by a factor inherent to 
the ASA-LMTO calculation, which we found to be under- 
estimated in ASA by 50%, compared to our full-potential 
result. Spherical approximations for the crystal poten- 
tial commonly underestimate the tendency to magnetism 
compared to the full-potential calculations in materials 
near a quantum critical point (cf., e.g., Ref. I32T ]. The 
other contributions calculated in LMTO are found via 
linear response formulas and do not include any renor- 
malization. Therefore, the final result in the rightmost 
column of the Table IIIII includes a renormalization by a 
factor of two for the LMTO core polarization, of 6.8 for 
all other LMTO contributions, and no renormalization 
for the LAPW results. 

We notice that the LAPW and the Stoner renormal- 
ized LMTO results for Ki SO (spin part) agree rather well; 
K lso = RKf s0 = -1.41 in LMTO and -1.44 in the non- 
relativistic LAPW (spin-orbit hardly affects this term). 
The small discrepancies can be attributed to small dif- 
ferences in the band structure and in the calculation 
procedure^ 

The dipole contribution to Ki SO is zero by symmetry. 
However, the orbital contribution is not, and it appears 
to be comparable to the contact term (as opposed to the 
magnetic susceptibility, where the orbital part is rather 
small). While nonrelativistic LMTO calculation cannot 
be used to evaluate this contribution, the WIEN2k code 
includes an option of relativistic calculations in an ap- 
plied magnetic field. Importantly, we found that it was 
absolutely necessary to include not only spin-orbit inter- 
action, but interaction of the orbital moment with the 
external field as well. We found that after including the 
orbital part, Ki SO — —2.56. Regarding the oxygen Knight 
shift, we get Ki SO = —0.09 and —0.05 for the in-plane and 
apical oxygen, respectively. 

Experimentally, Ki SO ss —2.98 was reported in Ref. |29|. 
and a slightly smaller value Ki SO w —2.91 in later papers 
by the same group 6 . Our calculated number of —2.56 is 
12% smaller than the later data, reminiscent of the 15% 
underestimation of the total susceptibility (Table HTj) . 36 
While the calculated Ki SO for oxygens has the wrong sign, 
this does not seem disturbing, since both the experimen- 



tal and the calculated numbers are very small and re- 
sult from a substantial cancellation of various anisotropic 
terms, and the absolute magnitude of the error is rather 
small. 

This agreement seems impressively good, but, just as 
with the magnetic susceptibility, this is probably decep- 
tive. As mentioned, Sr2Ru04 is close to a magnetic quan- 
tum critical point, and therefore belongs to a growing 
class of materials where LDA overestimates the tendency 
to magnetism (see, e.g., Ref. l37l and references therein). 
It is generally believed that magnetic fluctuations play 
an important part in linear response in such materials, 
suppressing the magnetic susceptibility compared to the 
essentially mean field LDA treatment 38 . Indirect evi- 
dence indicates that LDA does indeed overestimates the 
propensity to magnetism; for instance, in LDA Sr2Ru04 
is unstable against formation of a spin density wave at 
the nesting vector (1/3,1/3,0)^. It is hard to imagine 
that the same fluctuations that stabilize the paramag- 
netic state as opposed to the spin density wave ground 
state in LDA would not reduce Xi K (and also 1/T\T). 
Also, calculations predict Sr2Ru04 to be ferromagnetic 
at the surface^, while experimentally it is not the case^i. 
We conclude that the good agreement of the calculated 
X and Ki SO is largely due to a cancellation of errors: we 
overestimate the spin susceptibility, but underestimate 
the orbital susceptibility, as also evidenced by our failure 
to reproduce the sign of the anisotropy in \~ Note that 
underestimation of the orbital polarization, due to corre- 
lation effects, is very typical for d-shells in the transition 
metals. 

Let us now analyze the anisotropy of the Knight shift. 
Experimentally, Ishida et al^ reported for Ru K an j,, = 
(K z - K x )/3 = -0.23% and later -0.20% (Ref. 
Within the spherical approximation for the potential, 
the LMTO calculations (Table HlT|) give K anlso = -0.05 
(which is solely a dipole effect). On the other hand, a full 
potential approach (LAPW) results in K an i SO = —0.23 in 
the dipole part, which is reduced to —0.21 by the orbital 
effects. This excellent agreement reflects the fact that 
the lion share of the calculated anisotropy comes from 
the dipole part. Note that a possible overestimation of 
the contact interaction would not affect K an i S0 . 

It is instructive to analyze why LMTO underestimates 
the dipole contribution compared to the full-potential 
calculations, by as much as a factor of 4.5 (for Ru), de- 
spite the fact that it is is very accurate for the contact 
term. To this end, let us recall that for the dipole inter- 
action, the contributions from xy and yz/xz bands have 
opposite signs and tend to cancel each other. Keeping 
only diagonal terms, which dominate, we find 

jsdip 2 / r —3\ 

KZZ so =^*-^^[2N xy -(N vz + N xz )}. (5) 

Here N xy , N yz and N xz are orbital- projected density 
of states at the Fermi level (Table IIV[) per atom and 
(r~ 3 )i = J dr\4>i(r)\ 2 /r 3 , where <j>i(r) is the radial solu- 
tion of the Schrodinger equation with angular momentum 
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TABLE IV: Projected density of states (states/Ry/atom), obtained with the LMTO method (set up with three empty spheres^) 
and comparison with LAPW. In the last two columns |0 S (O)| 2 and (1/r 3 ); = f drr 2 \(f)i(r)\(l/r 3 ) (a.u.) are displayed, where 
4>i (r) is the radial solution of Schrodinger equation for the free atom and I its angular quantum number. 
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quantum number I. We found in ASA-LMTO (see Table 



W\) that N. 



N yz , and therefore the dipole con- 



tributions nearly cancel each other. On the other hand, 
this cancellation is far less complete in LAPW (see Ta- 
ble |IV|: in LMTO the ratio N xy /N xz w 1.2, while in 
LAPW it is nearly 2. Similar discrepancy is observed in 
the N x /N y /N z ratios for oxygens. On the other hand, 
the total DOSs, Nd and N p agree reasonably well. Al- 
though one can expect that the spherical approximation 
somewhat underestimates the DOS anisotropy, the ob- 
served discrepancies seem to be too large. A close in- 
spection reveals the secret: LMTO underestimates the 
electrostatic crystal field splitting, that is, the shift be- 
tween the on-site energies of the d xy and d xz orbitals, 
by approximately 4 mRy, which leads to a relative shift 
of the xy and xz/yz band by ±4 mRy compared to the 
LAPW calculations. This seems like a small effect, and 
it is, but because the Fermi level is situated on a slope 
of the DOS, shifting these bands in the opposite direc- 
tion leads to a serious redistribution of the DOS at Ep 
between them. The more accurate NMTO bands, with 
a better account for nonspherical effects, are closer to 
LAPW: we find N xy - N xz = 4.9 states/Ry/atom. 

It is more complicated to analyze K an i SO for oxygen, 
because of bigger experimental uncertainty, strong tem- 
perature dependence (on the 01 site), and somewhat 
larger orbital contribution in the calculations. The in- 
plane anisotropy on the oxygen bridging two Ru along 
x, defined as K\ niso = {K y - K x )/2, is -0.32% in 
the experiment, while in the calculation it is —0.21, 
mostly dipole. The out of plane anisotropy, K^ niso = 
[K x - (K x + K y )/2}/3, is 0.04 in the experiment and 0.01 
in the calculations. Again, one ought to keep in mind 
the overall small magnitude. Finally, for the apical oxy- 
gen K an i SO is essentially zero both in the experiment and 
theory. Note that in case of oxygen, LMTO also under- 
estimates the dipole contribution, by virtue of the same 
argument (for O, K z a -U%[2N Z - {N x + N y )] {r" 3 )^. 



Finally, strontium appears to be the only site where 



the Fermi contact term dominates. Unfortunately, to the 
best of our knowledge at present there are no experimen- 
tal data available for Sr. 



V. NMR RELAXATION RATES 



Let us now analyze the relaxation rates. We calculated 
1 JT\ with the procedure described in Ref. [33j and based 
on the LMTO band structure. Within this approach, 
the relaxation rate is obtained from Fermi's golden rule 
and linear response theory; for a polycrystalline sample 
1/T±T can be expressed as 



1 



Tr-|HiV| 2 



(6) 



where H is the hyperfine field operator previously intro- 
duced. Since the cross-terms in Eq. ^ vanish exactly for 
a polycrystal^, the relaxation rate is, apart from the core 
polarization contribution, the sum of contact, a dipole- 
dipole and an orbital term 

^ f = 2nk B h 1 ^ [|H c Af + |H d Af + |H°7V| 2 ] . (7) 

We calculated these contributions using the LMTO ba- 
sis set; the results are displayed in Table El The core 
polarization term was obtained separately from the core 
Knight shift and Korringa law. 

Available to us LAPW WIEN2k package currently 
does not allow for calculations of the relaxation rate. 
However, we can make use of the understanding gained 
on the previous stage by comparing LMTO and LAPW 
results for the Knight shifts to evaluate the reliability of 
the LMTO calculations. 

We obtained the bare relaxation rate (Table [Vj and 
calculated the Stoner enhancing factor separately. The 
computation of the Stoner enhancing factor for the re- 
laxation rate, Ri/^Ti requires knowledge of the q depen- 
dence of the spin susceptibility, according to the following 
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TABLE V: l/T x T in (K sec)" 1 , results based on LMTO calculations. The Stoner factors R 17 « 26 and 2 17 « 3.25 (core) are 
used for the renormalized values (see text). The 4th and 5th column display the total relaxation rate obtained after Stoner 
renormalization of, respectively, the spin terms and of both spin and orbital terms. 





orbital 


dipole 


contact 


core 


Tot (spin renorm) 


Tot (all renorm) 


Expt. ° 


Expt. b 


101 Ru 


0.18692 


0.02855 


0.00004 


0.53340 


2.67 


7.3 


15 




17 oi 


0.03291 


0.01392 


0.00007 


0.00541 


0.41 


1.23 


0.8 


1.1 


17 02 


0.00085 


0.00028 


0.00011 


0.00004 


0.011 


0.032 


0.025 


0.025 


87 Sr 


0.00029 


0.00009 


0.00254 


0.00003 


0.069 


0.076 







"For 101 Ru Ref. \MM (T=4.2 K). For 17 (planar and apical) 
data are from Ref.|4j(T=4.2K). 
frRef.El sample with T c = 1.5K (T = 4.2K). 



formula*^: 

\[l-7(q)i?exo(q)] / 

Using the Lindhart susceptibility we showed earlier— 
that in 3D, R\jT x T ~ R 5 ^ 3 - However, Sr2Ru04 is rather 
2D, in which case the real part of the noninteracting sus- 
ceptibility is constant, and, assuming a q-independent 
I, we get simply Rx/TxT ~ R 2 , as first pointed out by 
Shastry and Abraham4£ Actual dependence of I on q 
in Sr2RuC>4 is discussed in Ref. HD, where it was esti- 
mated as I(q) = I /[I + b(q/G) 2 ], where G = ir/a, and 
b is a numerical constant close to 0.08. Using this with 
Eq. ([SJ) and taking into account the actual q dependence 
of /mxo(q) we find Ri/^t ~ R 1 ' 7 = 26. We use this en- 
hancement factor Ri/^t to renormalize the contact and 
the dipole-dipole contributions to the relaxation rates. 
Similar to the Knight shift, as discussed in Section IV, 
the calculated core polarization contribution to the relax- 
ation rate is Stoner-renormalized according to the LMTO 
Stoner factor, which by itself is smaller than the renor- 
malization factor in LAPW. Thus, the core-polarization 
part was renormalized by a factor 2 17 . 

The renormalization of the orbital terms is more con- 
troversial. Even if we knew the renormalization factor 
for the spin shifts exactly, there is no guarantee that the 
renormalization for the orbital shifts is the same. Indeed, 
as discussed above in connection with the Knight shift, 
the orbital polarization is induced by both direct inter- 
action of the orbital currents with the external field, as 
well as, indirectly, via the spin-orbit interaction with the 
induced spin density. The former process is not a subject 
of Stoner renormalization, while the latter is. It is there- 
fore likely that the renormalization of the total orbital 
term is not as large as the Stoner renormalization of the 
spin terms. In Tab. |V]we show the total relaxation rates 
obtained by renormalizing the spin shifts only, and those 
obtained renormalizing (with the same enhancing factor 
w 26) the orbital shifts too. These numbers are the lower 
and upper limit for the LDA relaxation rates. Judging 
from the LAPW Knight shift (with spin-orbit) one may 
expect that the spin-orbit contribution to the relaxation 
rate is larger for for Ru than for O; correspondingly, the 



renormalization of the orbital terms is expected to be 
larger for Ru than for O. 

Since there are so many uncertainties in the renormal- 
ization of the orbital terms, first we discuss the relaxation 
rates obtained after renormalization of the spin (contact, 
core and dipole-dipole) terms only. The total calculated 
\jT\T » 2.67 (K sec) -1 , underestimating the low tem- 
perature \/T\T w 15 (K sec)" 1 by a factor 5. For 01 we 
find large orbital and dipole-dipole terms, and a sizable 
core polarization contribution. We find \ jT\T = 0.41 
(K sec)" 1 . Experimentally, \jT\T was found to be 0.8 
(K sec)" 1 at 4.2 K; now the calculated value is too small 
by a factor 2. For 02 the value l/T x T = 0.023 (K sec)" 1 
was measured4^i Our calculated value X/T\T = 0.011 
(K sec) -1 , again to small of a factor 2. 

In order to understand better the source of this dis- 
crepancy we examine the case of Ru in more detail. The 
leading contribution to 101 l/TiT is the core term; judg- 
ing from the Knight shift calculations, we do not expect 
the core polarization to be particularly poorly described 
by LMTO. The next largest term is the orbital contribu- 
tion. It is given, for Ru, for example, by the following 
expression 

GrV) b = \ C ^B{r-*)2f[Nl z + Nl 

+ 4N 2 xy + 2N xy (N xz + N yz )}, (9) 

where C = (4Trk B /h){~/ n /j e ) 2 w 0.404 x 10 5 /(Ks) for 
Ru 101 . Note that the orbital contribution to the re- 
laxation rate does not require spin-orbit coupling (it ap- 
pears already in the nonrelativistic approximation), but 
the spin-orbit coupling can and does contribute to the 
orbital polarization. Note that if the spin-orbit contribu- 
tion to the relaxation rate would be of the same relative 
size as the corresponding contribution to the Knight shift 
(that is, of the order of the core contribution) we would 
expect an additional spin-orbit induced relaxation rate of 
about w 1.75 (K sec) -1 . This would be the first contri- 
bution missing in Table fVl 

Furthermore, the orbital term is also partially under- 
estimated because of the difference between the LMTO 
and LAPW DOSs, as already discussed in the previous 
session. Using LMTO partial DOSs in the approximate 
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formula Eq. © we find {l/TiT) orb w 0.29 (K sec)" 1 , 
in good agreement with the full calculations (Tab. W\ : 
using instead LAPW DOSs (the exact comparison of 
LAPW and LMTO is of course not possible because space 
is divided up differently in the two methods) we find 
(l/TiT)orb w 0.4 (K sec) -1 , which indicates that the 
LMTO orbital term is understimated by a factor w 1.4. 
Similar considerations apply to the dipole-dipole term, 
described by 

+ W* y +3N xy (N xz + N yz )]. (10) 

This correction (for both terms) is about 0.4 (K sec) -1 . 

Both these corrections together raise the estimate for 
the relaxation rate on Ru to ~ 5 (K sec) -1 , still short 
of the experimental number of 7.3 (K sec) -1 . We take it 
to be yet another manifestation that the orbial effects on 
Ru are underestimated in LDA and that the the leading 
relaxation mechanism is orbital. 

One should expect that correlation effects on O are 
rather moderate and then, after correcting for the 
LAPW-LMTO DOS differences, we should get a reason- 
able agreement with the experiment (spin-orbit effects 
are also weak for O). Indeed, after this correction we ob- 
tain (l/TiT) i « 0.65 (K sec) -1 and (l/TiT) 02 w 0.16 
(K sec) -1 , reasonably close to the experiment. 

Finally, we neglected the quadrupolar contribution to 
1/TxT. For 101 Ru this term is small, but not negligible^; 
taking this term into account would therefore increase 
the agreement with experimental measurements. Inter- 
estingly, the quadrupolar term could partially contribute 
to smear out the Hebel-Slichter peak. The quadrupo- 
lar contribution is however negligible for 01 and 99 Ru; 
measurements of the relaxation rate below T c for these 
two ions could therefore shade further light on our un- 
derstanding of NMR data and clarify the absence of the 
Hebel-Slichter peak. 

VI. SUPERCONDUCTIVITY 

Finally, we want to make a few comments on pos- 
sible alternative interpretation of the invariance of the 
Knight shift across the superconducting transition. As 
we discussed in the introduction (see also discussions in 
Refs. HUTU) the interpretation in terms of the order pa- 
rameter, rotating in the field of 0.02 T, less that one per 
cent of the superconducting gap, is very unlikely. The 
interpretation of the temperature independence of the 
Knight shift in the superconducting state of Sr2Ru04 
in terms of a chiral p-state is not possible, at least not 
on the level of existing theory. Given that this is the 
only experiment uniquely identifying the superconduc- 
tivity in Sr2Ru04 as chiral p-wave, this opens the door 
for other possibilities regarding the symmetry of super- 
conducting state, such as the chiral d-state proposed in 



Ref. 1, A oc xz i iyz, or a planar p-wave state, as sug- 
gested in Ref. [TH 

Given the rather inexplicable absence of the Knight 
shift decay below T c , it is tempting to bring up an anal- 
ogy with such materials as V and Hg. Two explana- 
tions of the same phenomenon in these elemental met- 
als have been proposed. The first one is the spin-flip 
scatterin g 16 ) 17 by the grain boundaries or point defects 
due to the spin-orbit coupling. Spin-orbit indiced split- 
ting of the conductivity bands is rather nonuniform over 
the Fermi surface, reaching 200K at some places (Figs.[U 
Anderson^ 7 - has derived a formula describing this ef- 
fect: the ratio of the Knight shift at T = and at T > T c 
is K sup / K norrn 7 s .j./6A, where 7 S ./. is the rate of the 
spin-flip events and A is the superconducting order pa- 
rameter, and 7s./. > A. If "f s .f. gets noticeably larger 
than A, the Knight shift changes rather little across T c . 
In principle, there is nothing unimaginable in "f s .f. > A, 
however, we should remember that nonmagnetic impu- 
rities in Sr2Ru04 are pair breakers, and samples with 
T c « 1.4 K, and the experiments we strive to explain 
were performed on such samples, have to satisfy the pair 
breaking condition j tr < A, where jtr is the transport 
relaxation rate. Even given the sizeable spin-orbit in- 
teraction in Sr2Ru04, it is totally unrealistic to assume 
spin-flip relaxation to be larger than momentum relax- 
ation, that is 7s./. <C 7tr < A should hold. 

Another mechanism that can emulate a temper- 
ature independent Knight shift was described by 
MacLaughlinA 9 : if the Fermi contact and the core po- 
larization contribution cancel each other, the remaining 
orbital part, being an entirely charge effect, will stay un- 
supressed below T c . This was supposed to be operative, 
for instance, in V. There are two objections against ap- 
plying the same scenario to Sr2Ru04. First, in Ru we 
find that the Fermi contact interaction is negligible, but 
the core polarization is responsible for about a half of the 
total Knight shift. Even though it is probably overesti- 
mated, we can guarantee it is not negligible. Moreover, 
for both 01 and Ru we find a sizeable dipole term (in 
fact, for 01 it is the leading term). In Ref. [l^ this term 
was neglected, but it also originates from the spins of 
electrons and should be sensitive to pairing. So, on this 
(probably oversimplified, as discussed below) level nei- 
ther of the two known explanations of the invariance of 
the Knight shift below T c is applicable. 

While this simplistic picture^ does not appear to work 
for Sr2Ru04 it brings up an important point, that the 
standard theory of the Knight shift in the superconduct- 
ing state, developed for Fermi-contact interaction in one- 
band superconductors has to be revised to be applicable 
to Sr2Ru04. Indeed, while the orbital interaction does 
not dominate the Knight shift on Ru, it constitutes at 
least a half of the total K, and probably as much as 
2/3, if we take into account that the core polarization 
is probably overestimated in LDA. This alone raises our 
expectation for K(T = 0) from zero to about 2/3 of its 
normal value. On top of that we have a sizeable dipole 
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term, which for the in-plane field is opposite in sign to the 
core polarization. Most interestingly, the temperature 
dependence below T c should be highly nontrivial even 
for a singlet pairing: as discussed, this term comes about 
through near-cancellation of the contributions from the 
xy and xz + yz bands. For several reasons, these bands 
are expected to react differently to magnetic field (and 
the actual measurements were performed in a sizeable 
field): they have different Fermi velocities along the c- 
axis, and possibly rather different gap values. Therefore, 
the partial densities of states decay (in a singlet case) 
below T c at different rates, with a possibility of the net 
dipolc Knight shift to decrease, increase or (accidentally) 
stay constant! The exact answer requires a quantitative 
analysis that goes beyond the framework of this paper, 
but we want to emphasize that the fact that the dipole 
contribution is strong calls for an entirely new theory of 
the Knight shift below T c . A similar consideration applies 
to the Knight shift on O. 

VII. CONCLUSIONS 

We have calculated from first principles the magnetic 
susceptibilty, the Knight shifts and the relaxation rates in 
Sr 2 Ru0 4 . We adopted different LDA techniques (LMTO 
and LAPW) and performed calculations with and with- 
out explicit account for the spin-orbit interaction. The 
agreement with available experimental data is good, in 
some cases surprisingly good. The isotropic magnetic 
susceptibility, the isotropic and anisotropic Knight shift 
for Ru are all reproduced better than within 15%. How- 
ever, there are notable discrepancies, of which most im- 
portant is, probably, the incorrect anisotropy of magnetic 
susceptibility. We argue that this error is related to an 
enhancement of the orbital polarization due to correla- 
tion effects beyond LDA. Despite the good quantitative 



agreement with the experiment of the isotropic suscpti- 
bility and isotropic Knight shift on Ru, we do think that 
we witness a cancellation of two errors: On one hand, 
we overestimate the spin susceptibility, which should be 
reduced by spin fluctuations missing in the LDA, and on 
the other hand we underestimate the orbital polarization 
in the Ru d shell. If this interpretation is correct, one 
has to conclude that orbital effects are extremely impor- 
tant and probably dominate both the Knight shift and 
the NMR relaxation on Ru. 

We have considered the accepted mechanisms of the 
Knight shift invariance in elemental metals and con- 
clude that they cannot explain the same phenomenon 
in Sr 2 RuC>4. In particular, we can exclude the spin-flip 
scattering by defects in Sr2Ru04. On the other hand, 
orbital polarization of electrons at the Fermi level, which 
was used earlier to explain a similar effect in vanadium, 
is also the leading contribution in Sr 2 Ru04. Although 
it does not entirely dominate the Knight shift, as in V, 
it does substantially reduce the expected effect of singlet 
superconductivity on it. Finally, sizeable contribution of 
the dipolc-dipole interactions together with the orbital- 
dependent order parameter essentially renders the exist- 
ing theory of the Knight shift suppression by a singlet 
pairing incomplete and calls for its revision. We hope 
that this conclusion will inspire further experimental and 
theoretical studies. 
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